NASA Technical Memorandum 100203 


A Detailed Description of the Uncertainty 
Analysis for High Area Ratio Rocket 
Nozzle Tests at the NASA 
Lewis Research Center 


Kenneth J. Davidian 
Lewis Research Center 
Cleveland, Ohio 

Ronald H. Dieck 

Pratt & Whitney 

West Palm Beach, Florida 

and 

Isaac Chuang 
Lewis Research Center 
Cleveland, Ohio 


Prepared for the 

24th JANNAF Combustion Meeting 
Monterey, California, October 5-9, 1987 


NASA 

<NASA-TH- 100203) A DETAILED DESCRIPTION OF N87-2E602 

TEE 0NCERTAINTY ANALYSIS FC£ klGB AREA RATIO 
POCKET KCZZLE TESTS AT THE NASA IEH1S 

RESEARCH CENTER (NASA) 30 p Avail: NIIS Onclas 

EC AC3/EF A0 1 CSCL 21H G3/20 00S7E53 


A Detailed Description of the Uncertainty 
Analysis for High Area Ratio Rocket Nozzle Tests 
at the NASA Lewis Research Center 

Kenneth J. Davidian Ronald H. Dieck 

NASA Lewis Research Center Pratt & Whitney 

Cleveland, Ohio West Palm Beach, Florida 

Isaac Chuang* 

NASA Lewis Research Center 
Cleveland, Ohio 


ABSTRACT 

A preliminary uncertainty analysis has been performed for the High Area Ratio Rocket Nozzle 
test program which took place at the altitude test capsule of the Rocket Engine Test Faciltiy at the 
NASA Lewis Research Center. Results from the study establish the uncertainty of measured and 
calculated parameters required for the calculation of rocket engine specific impulse. A generalized 
description of the uncertainty methodology employed is provided. Specific equations used and a 
detailed description of the analysis are presented. Verification of the uncertainty analysis model was 
performed by comparison with results from the experimental program’s data reduction code. Final 
results include an uncertainty for specific impulse of 1.30%. Largest contributors to this uncertainty 
were calibration errors from the test capsule pressure and thrust measurement devices. 

INTRODUCTION 

The specific impulse of a rocket engine is one of the most important performance parameters. 
Results of mission analyses and vehicle sizing studies are dependent on the accuracy of the value 
used for an engine’s specific impulse. Computer predictions of specific impulse for engines with 
high area ratio nozzles (a domain for which the codes were not designed) can be compared with 
experimentally measured values to verify the analytical output. To insure meaningful correlation 
between computer and test results, the specific impulse should be measured as accurately as possible. 
This report details an uncertainty analysis performed on data from rocket engine tests at the NASA 
Lewis Research Center. 

A measurement uncertainty methodology for chemical rocket engines has been establish by JAN- 
NAF (Joint Army Navy NASA Air Force), formerly the Interagency Chemical Rocket Propulsion 
Group (ICRPG) (ref. 1). A goal of 0.25% on specific impulse uncertainty has been suggested and 
has been the topic of recent study (ref. 2). 

The objective of this study was to determine the uncertainty on the calculated value of specific 
impulse from measurements taken during the High Area Ratio Nozzle test program. Uncertainties 
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for measured parameters used in the calculation of specific impulse are given along with thoses for 
other measured and calculated parameters of interest. 

The following assumptions were made in the analysis: 

1. The experiment for which the uncertainty analysis was performed consists of carefully con- 
trolled measurement processes and calibration corrections for all instrumentation were per- 
fectly carried out. This eliminated large bias errors of known magnitude. Small bias errors 
were assumed to be zero. In reality, the bias errors of measuring devices are not zero, although 
they generally are very small. An example of a small bias of unknown magnitude would be the 
reported uncertainty of a transducer which was calibrated by the National Bureau of Standards 
and which was incorporated into the test set-up. The level of uncertainty is unknown (only 
the range is known) and is fixed. Hence, it contributes to the experimental uncertainty as a 
bias error. 

2. All data reduction errors were zero. For this experiment, this assumption is reasonable since 
the computer truncation errors and high order curve fit errors are negligible. 

Although a general expectation of uncertainty may be obtained from an uncertainty analysis, 
many uncertainty analyses are tailored specifically for a given test, a given set of instruments, or 
a given set of data reduction equations. For this study, the experiment being analyzed is a rocket 
engine which employs gaseous hydrogen and gaseous oxygen as its fuel and oxidizer. The experiment 
was run at the NASA Lewis Research Center’s Rocket Engine Test Facility (RETF) inside a low 
pressure capsule to simulate altitude conditions. 

GENERALIZED DESCRIPTION OF THE UNCERTAINTY METHODOLOGY 

The uncertainty methodology is described for different applications in references 1 and 3 through 
27 and is described in detail in references 3 and 19. Only the basic principles will be discussed here 
for clarity. 

Measurement Error 


All measurements have uncertainties which are associated with experimental errors. These are 
the differences between the measured and an accepted standard true value. The total error of a 
measurement is usually expressed in terms of two components: a random (precision) error and a 
fixed or systematic (bias) error. 

Precision Error 

The distribution of a set of data points caused by random error are characterized by the standard 
deviation. The precision error range is defined to be twice the standard deviation for a large sample of 
data (more than 30 data points) and includes approximately 95% of the total scatter of measurements 
(ref. 27). The precision index of the data, S*, approximates the true value of standard deviation 
and is defined by: 

% = (*> 

where N is the number of data points averaged, A' is the average of the N data points, and A\ is 
the i th individual data point. 
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The precision index of the average of a set of measurements is always less than that of an 
individual measurement. The precision index of an average, 5^, is defined as: 


c. l / eL (^-^) 2 

x ~V N(N-l) 


( 2 ) 


It should be noted that the precision index, Sx , may also be estimated by assessing the variability 
of data about a curve fit. This variability, called the Standard Estimate of Error (SEE) is calculated 


as: 


SEE = 


N -C 


( 3 ) 


where y x is the i th y value, y iiC is the equivalent y t calculated from the curve fit at X x} N is the 
number of data points, and C is the number of constants in the curve fit. 

If many estimates of the same precision index are available, a better estimate of Sx may be had 
by “pooling” the estimates as follows: 




E* 1 *5? 

N 


( 4 ) 


where S p is the pooled precision index and v x is the degrees of freedom of the i th estimate of S x . 

Bias Error 


Bias error is a systematic error which remains constant during a given test. Thus, in repeated 
measurements of a given variable, each measurement has the same bias. Note that it is assumed 
that the uncertainty analysis is applied to a carefully controlled measurement process within which 
all known calibration corrections have been made. If these calibrations were perfect, there would be 
no bias error. However, there always remains some systematic error of unknown magnitude. This 
bias error remaining must be estimated. If there is no statistical equation with which to determine a 
value, the estimate of the bias limit, B, must be based on judgement. Reference 27 provides several 
methods for estimating this bias error. 


Combining Errors 

In describing the precision and accuracy of a measurement, there is a need for a single value to 
quantitatively categorize the data adequacy. In obtaining this value, it is useful to group the error 
source. 

Sources of error can be divided into three categories: calibration errors, data acquisition errors, 
and data reduction errors. Each of these sources of error have components of bias and precision 
error. 

To obtain the bias of a given parameter (such as temperature or pressure), the root sum square 
(RSS) method is used to combine the bias limits from numerous ( N ) elemental sources of error. 
Thus 

b = Jb[+bi + --- + b% (5) 

where B is the bias limit of a parameter and Bi is the bias limit for the i th elemental error source. 
Similarly, 

5 = ^/s? + Si + ■■■+ S 2 N (6) 
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where S is the precision index of a parameter and 5, is the precision index for the i th elemental 
error source. 

For most experimental decisions, a single number is needed to express a reasonable limit of error 
for a given parameter. A model for combining the bias and precision errors, therefore, must be 
adopted which will yield the interval X ±U i where U is the measurement uncertainty. Z7 is a band 
within which the true value of the parameter is expected to lie, for some specified coverage. While 
no rigorous confidence level can be associated with the uncertainty, coverages analogous to the 95% 
and 99% confidence levels can be given for the two recommended uncertainty models (ref. 27). Thus 
one uses 

r/99 = B + t^S (7) 

for approximately 99% coverage and 

u 9 5 = v/B 2 + (< 95 S ) 2 ( 8 ) 


for approximately 95% coverage. { 7 99 and C/ 9 5 are also referred to as Uadd and Urss, respectively. 

The Student’s t value , t 95 , is a function of the degrees of freedom used in calculating 5*. For 
large samples, (i.e., N > 30), <95 is set equal to 2. Table I gives the values for <95 for different values 
of degrees of freedom. Otherwise the Welch-Satterthwaite formula is used to provide the degrees of 
freedom, v, according to 


. = <z£itf )! 

2^t=i Ui 


( 9 ) 


where S* represents the precision indices of the various error sources involved, and v t represents the 
degrees of freedom associated with those error sources. 


Uncertainty o f a Result 


Errors in measurements of various parameters (P) are often propagated into a derived result 
(r) through the functional relationship between the result and those parameters. The existence of 
such a relationship requires the use of influence coefficients (0* ) which are used to obtain the error 
propagated to the result because of a unit error in the parameter. Thus if 


r = /(/>!, P 2 ,P 3 ,...,P W ) 


(10) 


where N is the number of parameters involved, and r is the computed result, then the influence 
coefficients are: 



(ii) 


It is important to note that, as in all uncertainty analyses, the bias and precision errors of the 
parameters are kept separate until the last step, computing the uncertainty of a result. 

The precision index of a result, S r , is given by 


JV 


Sr = 


T>s *) 2 


1=1 


Similarly, the bias limit of a result, B r , is given by 

B r = 


1 


N 




i = l 


( 12 ) 


(13) 
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The uncertainty of a result is given by either of the two uncertainty models shown above. The 
Student’s t value is defined as before. 

The Welch-Satterthwaite formula is used to provide u T and is slightly more complex. It’s formula 


is 


f ) 4 


(14) 


Note the influence coefficients convert each Sp i into units identical to those of the result, a necessary 
step for combination. 


UNCERTAINTY ANALYSIS DESCRIPTION 


Facility Description 

The altitude test chamber in the RETF includes a test capsule, diffuser, spray cooler, ejectors, 
liquid drain lines, and the water detention tank. The exhaust gases of the rocket engine aid in 
altitude pumping by passing through a second throat diffuser before exhausting into the spray 
cooler. Approximately half of the exhaust gases are condensed to a liquid and pass down the drains 
to the water detention tank. Ejectors, driven by gaseous nitrogen, pump the remaining exhaust 
gases through two short stacks to the atmosphere. 

The engine instrumentation system is displayed in Fig. 1. Propellant flow rates are determined 
using calibrated venturis. Temperatures are measured using Chromel-Constantan thermocouples. 
A thermocouple type vacuum gauge is used to measure the vacuum reference pressure while the 
remaining pressures are measured by strain-gauge bridge type pressure transducers. Absolute and 
differential pressure transducers are used. The thrust stand is capable of measuring thrust levels to 
13.3 KN (3000 Ibf) and was designed to have a random error of less than ± 0.1% of full scale. For 
more information on the facility, see reference 28. 

Instrumentation in the facility’s data acquisition system provides analog signals that are recorded 
and converted to a digital signal by an automatic data digitizer at a rate of 50 readings per second 
per parameter and sent to an IBM 370 computer. The computer averages the values in groups of 
five to provide data output at ^ second intervals. 

Site thrust, the force actually acting on the test stand, was measured by a triad of load cells in 
a parallel configuration to account for a thrust vector which is not perfectly aligned with the rocket 
engine’s centerline axis. Similarly, other transducers may be shown in Fig. 1 as one device when, 
in fact, there may be more than one instrument recording data at that location. In these cases, the 
redundancy was to identify faulty transducers should one fail. 

Upstream of the subsonic venturi, propellant line pressure and temperature were recorded. Just 
upstream of the injector, each propellant’s pressure and temperature was again measured. Pressure 
transducers were located in the combustion chamber and differential transducers between the in- 
jector and chamber pressure transducers to allow for direct measurement of injector pressure drop. 
Redundant measurements of injector pressure drop were also obtained by subtracting the chamber 
pressure from the injection pressure. Differential pressure transducers, referenced to a vacuum tank, 
were used to measure nozzle wall static pressure and capsule (ambient) pressure. 

Equations Used to Reduce Data 


Experimentally, only the measurements of site thrust, capsule pressure, propellant line pressure 
and temperature, and differential pressure from the propellant line to the venturi throat were required 


5 



to calculate the rocket engine’s specific impulse. Data reduction equations used to compute specific 
impulse are described below. 

Pressu re at the ven turi throat. Throat pressure of the hydrogen and oxygen subsonic venturi 
flow meters, pa V) was calculated by: 

Pth = Pixne - dp (15) 

where pu ne is the static pressure measurement upstream of the venturi and dp is the differential 
pressure measurement from the pane transducer to the venturi throat. 

Line contract i on ratio . Propellant line contraction ratio, /?, was defined by: 


0 = 


dth 


dline 


(16) 


where d t h is the venturi throat diameter and di xne is the propellant line diameter. 

Thermodynamic properties. Using the measured values of fluid pressure and temperature up- 
stream of the venturi throat, the enthalpy, density, and entropy of the propellants in the propellant 
lines were determined using the Gas Properties program (GASP) (ref. 29). By assuming isentropic 
flow to the venturi throat, the throat pressure, calculated using Eq. (15), and the propellant entropy 
values were input to GASP to determine the fluid temperature, density, and enthalpy at the venturi 
throat. 

Velocity at the venturi throat. Velocity of each fluid at the venturi throat, V t k , was determined 


by: 


V th = 


f K (hj, ne h>th) 


(17) 


where K was a conversion constant equal to 50079.5 [ btu/ik L ^Une and h t h are the propellant line 
and venturi throat enthalpies, and pu nt and p t h are the propellant line and venturi throat densities. 
Propellant mass flow rate . Mass flow rates for each propellant, m, were determined from: 


m= C d jd 2 th p th V th (18) 

4 

where Cd is the venturi discharge coefficient. Discharge coefficient error was determined by inde- 
pendant calibration and its calibration incorporates any error in the determination of venturi throat 
diameter. 

The set of equations presented to this point, with the exception of Eq. (16), were employed four 
times in the data reduction procedure, once for each of the two oxygen line differential pressure mea- 
surements and once for each of the two hydrogen line differential pressure measurements. Hence, two 
values of oxygen and hydrogen throat pressure, thermodynamic properties, venturi throat velocity, 
and mass flow rates were calculated. Both oxygen mass flow rates were combined to arrive at an 
average value of oxygen mass flow rate, and the same procedure was followed for the fuel mass flow 
rates. 

Total mass flow rate . The total mass flow rate, rh TC > T , was simply the sum of the mass flow 
rates of each fluid, namely: 

m TOT - m ox + m fu (19) 

Vacuum thrust . Since the site thrust, F,iu, was measured in an imperfect vacuum, a correction 
was applied: 

Fvac = F tilt “1“ Pamb^-exit (20) 

where F vac is the vacuum thrust, p am b is the ambient, or capsule, pressure and A exit is the rocket 
nozzle exit area. 
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Vacuum specific impulse , 
rameter was calculated by: 


The performance parameter, vacuum specific impulse, I tp . 


I* P — 


F 

* vac 

TTlroT 


This pa- 

( 21 ) 


Propagation of Errors 


Errors which occur in measured parameters are propagated to the calculated parameters through 
the use of influence coefficients. Sixteen basic measurements, or primary parameters, were collected 
in this experiment. These parameters were the propellant line pressures, temperatures, and diame- 
ters, differential pressure from the line to the subsonic venturi throat, the venturi throat diameter, 
the nozzle exit area, the measured thrust, and the altitude capsule (ambient) pressure. 

These quantities were used to calculate certain variables, secondary parameters, such as the 
geometric and thermodynamic quantities. In turn, these were used to calculate the tertiary, or 
performance, parameters. Influence coefficients for the tertiary parameters with respect to the 
secondary parameters were calculated. However, these were not used to propagate the errors from 
the secondary to the tertiary results. Chain rule applications were employed to arrive at the partial 
derivative of the tertiary parameter with respect to the primary parameters. Since the secondary 
parameters were often functions of the same primary parameters, the partial derivatives of tertiary 
to secondary parameters did not accurately describe the ultimate effects of primary parameters upon 
the tertiary quantities. 

To validate the values of the partial derivatives which were used for the uncertainty analysis, 
inputs to the data reduction program used for the rocket engine tests were parametrically perturbed 
to determine the effect this would have on calculated specific impulse. The result was a value for the 
change of specific impulse with respect to a change of one of the primary parameters. By comparing 
output from the data reduction computer program with the partial derivates computed using the 
uncertainty analysis model, the validity of the influence coefficients, as well as the uncertainty 
model as a whole, was confirmed. An attempt was made to use the partial derivatives of tertiary to 
secondary parameters and calculate the effect of primary variables upon tertiary results through the 
error propagation equation, Eq. (12). Values for the influence coefficients did not match the data 
reduction program results, indicating the incorrectness of the attempt to calculate errors of tertiary 
parameters using intermediate (secondary) results. This mistake was basis for the author’s previous 
work (ref. 30) but has been corrected for in this study. 

For example, the propellant line density and enthalpy are calculated using the line temperature 
and pressure. Errors are propagated from the temperature and pressure measurements to these 
thermodynamic properties using influence coefficients. Later, the venturi throat velocity is calcu- 
lated as a function of the propellant line enthalpy and density. Velocity error is a combination (in 
part) of the propellant line enthalpy and density errors, however, since both of these properties are 
functions of the same basic measurements (temperature and pressure), combination of errors at the 
thermodynamic properties level leads to a greatly inflated value of error for the venturi throat ve- 
locity value. Therefore, the influence of propellant line pressure and temperature on velocity at the 
venturi throat, as described by partial derivatives calculated using the chain rule, must be computed 
to insure that all variables are independent of each other. 

Description of the Influence Coefficients 


Influence coefficients are required in the uncertainty methodology to propagate errors in measured 
parameters to errors in calculated parameters. Wherever possible, the influence coefficients were 
derived by using partial differentiation on the data reduction equations. To calculate the influence 
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coefficients for the propellant line and venturi throat thermodynamic properties, a perturbation 
method was used. Where partial differentiation of a calculated parameter was taken with respect to 
another calculated parameter, the chain rule was applied to attain the partial derivative with respect 
to one of the measured parameters, i.e. line pressure, line temperature, and differential pressure. 
Calculations for all influence coefficients are described below. 

Pressure at the venturi throat . Computation of these influence coefficients was accomplished 
easily: 

dpth 


dpiine 

tyth 

ddp 


= 1 


= -1 


( 22 ) 


(23) 


Absolute values of these coefficients are unity for both propellants. This results from the fact that 
there exists a one-to-one relationship between an error in line pressure or differential pressure to an 
error in calculated throat pressure. Units of these influence coefficients are psia/psia. 

Line contraction ratio . The line contraction ratio influence coefficients were found to be: 


d/3 _ 1 

dd-th di tne 


(24) 


d(3 _ dth 

ddline df ine 


(25) 


Although the propellant line diameters were the same for both propellants, the oxygen and hydrogen 
venturi throat diameters were different which resulted in slightly different values for the influence 
coefficient expressed by Eq. (25). Typical values are 0.62 inches’’ 1 for the value of Eq. (24) (both 
oxygen and hydrogen terms) and -0.16 inches _1 and -0.17 inches -1 for the value of Eq. (25), for 
the oxidizer and fuel lines, respectively. 

Thermodynamic properties. Influence coefficients of each calculated property with respect to 
both input properties were computed by perturbing each input parameter by a small amount (typ- 
ically 1%) and noting the effect that this change produced in the result. Tables II and III show 
typical values for each influence coefficient calculated in this manner. 

Since the partial derivatives in Table III are taken with respect to calculated parameters, they 
need to be reduced via the chain rule to partial derivatives with respect to line pressure, line temper- 
ature, and differential pressure. Equations (26), (27), and (28) give the generalized expressions which 
express the errors in the venturi throat thermodynamic properties as a function of the measured 
quantities. 


d</> 

dpitne 


d<f> dpth d<t> doling 
dpth Spline d&lint ^Pltne 

d<t> _ d<t> &Pth 
ddp dp th ddp 


(26) 

(27) 


&4> _ &<f> d*ixnt 

dTune dsu ne dTu ne 


(28) 


where <f > , the generalized variable, represents venturi throat enthalpy, density, and temperature. 

Velocity at the venturi throat . Expressions for the partial derivatives of the venturi throat pro- 
pellant velocity, Eq. (17), are shown below. 


_avth_ Jk_ 

dhl ' n ' 


8 



( 30 ) 


dVth = BV tH 
dhth dhline 
dVth _ __pj^ VtHp 

dpiinc pf ine [1 - 

dVth _ Pth VthP 

dVth (pth _\ 2 Wth£ 

d/3 U«J 


(31) 

(32) 

(33) 


Typical oxygen an 


-j i — i i:«« — r~.. t?~ (oo) i no /*/ <ec 

.AlU, UJUIV^CII lliic TOIUW IU 1 \AJSSJ tMV xv« BTUjlb 


o n /J IQ 

BTUflb „ 


/*/* 


respec- 
tively. A factor of five difference between the propellant lines’ values results from the large enthalpy 
decrease in the fuel line as compared to the oxygen line. Typical values of Eq. (30) are simply the 
negative of those for Eq. (29). Eq. (31), the partial derivative of throat velocity with respect to 
line density, has typical values of -0.24 and -"31 for the oxygen and hydrogen lines, 

respectively. The two orders of magnitude difference results from two causes. One is the throat 
velocity term in the numerator of Eq. (31). Oxidizer velocity at the venturi throat is 250 ft/ sec 
compared to the fuel line’s 1340 ft/sec } a factor of five difference. The other cause is that the 
density expression in Eq. (31) has a value of 0.24 for the oxygen line and 6.85 for the hydrogen 
line, approximately a factor of 20 difference. Typical values for Eq. (32), the partial derivative of 
throat velocity with respect to throat density, are different by two orders of magnitude also (0.25 


ib^ /ft* * OT ox y8 en l ine an d 32 for hydrogen line) for the same reasons mentioned above. 

Lastly, the throat velocity partial derivative with respect to the line contraction coefficient, Eq. (33), 
has typical values of 7.5 and 49.5 for the oxidizer and fuel lines, respectively with units of ft/ sec. 
The factor of five difference is primarily a result of the comparison of the values of throat velocity 
between the two propellant lines. 

To attain the partial derivatives of venturi throat velocity with respect to the measured quantities, 
the following equations are applied: 


dVtk = dv th dh hne | dv th dh th | dv th d Plxnt | dv ih d Pth 

dpitne dhltne &Pltne dh th &Pl\ne dpltne &Plint dpth &Plint 

dVth _ dv th dhth BVth Bpth 

ddp ~~ dh th ddp dpth ddp 

dVth _ dVth dhjjne dV t h dh th dV t h dpu nt dV t h dp th 

BTiine dhiine dTiine dhth dTiine dpHne dTi ine dpth dTn ne 

dV th = dVth d (3 

ddn ne d/3 ddimc 


(34) 

(35) 

(36) 

(37) 


dVth _ dVthd(3_ 

ddth d/3 dd th K } 

Values for the oxygen and hydrogen line partial with respect to line pressure were -0.19 
and -1.05 Equation (35) results in values of 4.9 and 16.5 for the oxidizer and fuel lines, 

respectively with units of . The partial derivative with respect to line temperature yielded 

values of 0.28 ^ e ^ - c and 1.25 for the oxygen and hydrogen lines. Equation (37), the line 
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diameter partial gave oxygen and hydrogen values of -1.2 and -8.3, and Eq. (38), the venturi 
throat diameter partial, resulted in values of 4.6 and 30.7 for the oxidizer and fuel propellant lines, 
respectively, all having units of . 

Propellant mass flow rate . There are three partial derivatives of the mass flow rate of each 
propellant. 

TT7 ~ + ‘M'thVth) ( 39 ) 

oa t h 4 Oath 

= + («) 

&-*!*** ,4,) 

Eq. (39), the partial derivative of propellant mass flow rate with respect to the venturi throat 
diameter, has typical values of 4.0 **^7/ and 1.3 for the oxygen and hydrogen lines, respec- 

tively. A factor of five difference exists between the oxygen and hydrogen line values of Eq. (40), the 
partial derivative with respect to throat density (the values are .21 and 1.4, respectively with units 
of due to the throat velocity term in the numerator of the expression. The last propellant 

mass flow rate partials is the throat velocity influence on propellant mass flow rate. Typical values 
are 0.003 for the oxygen line and 0.0002 for the hydrogen line with units of **7/17/ ■ The ver y l ar 8 e 
value of the oxygen line density is the cause of the difference between values. 

As in the case of venturi throat velocity, the influence coefficients of the propellant mass flow rate 
need to be expressed with respect to the quantities which were actually measured. The following 
expressions were used to calculate the required coefficients: 


dm __ dm dV th dm dp th 

dpiine dV th dpu ne dp t h dp\ ine 


(42) 


dm _ dm dV t h dm dp t h 
dTi ine dV t h dT line + dpth dTi ine 

dm _ dm dV t h dm dpth 

ddp dV th ddp + dpth ddp 
dm _ dm dV t h 
ddn ne dV t h ddn ne 


(43) 

(44) 

(45) 


Values for the line pressure partial, Eq. (42), were 0.0007 and 0.0002 for the oxidizer 

and fuel lines, respectively. Equation (43), the partial with respect to line temperature, had values 
of -0.0009 and -0.0002 for the oxygen and hydrogen propellants, respectively. The partial 

derivatives with respect to differential pressure took values of 0.015 and 0.0032 * or ^e ox yg en 

and fuel side, respectively. Each propellant's mass flow rate partial derivative with respect to 
propellant line diameter have values of -0.0039 and -0.0018 ■ 

Total mass flow rate . Like the throat pressure equation, the partials for the total mass flow rate 
were simply calculated as: 


dm TOT 

dm ox 

dmTOT 

dm fu 


= 1 

(46) 

= 1 

(47) 
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This results from the purely additive nature of the calculation, Eq. (19). Calculating the influence 
coefficients which relate back to the propellant line temperatures, line pressures, and differential 
pressures result in: 


dmpQp draper dra 

dpiint dm dpiinc 

(48) 

dthpor draper dra 

8T Unt dm dTu nt 

(49) 

draper drapeT dra 

ddp dm ddp 

(50) 

Am rim Am 


~ J Ui _ ^ ''l w X v - 

ddi xnt dm ddii nc 

(51) 

draper dmper 

(52) 

dd t h dm dd t h 

Each of these equations must be applied to the oxygen and fuel lines. 

Since the values of Eqs. 

(46) and (47) are unity, the values for the expressions are the same as those given for Eq. (39) and 

Eqs. (42) through (45). 

Vacuum thrust. Influence coefficients for the vacuum thrust equation 
and are expressed below. 

were simple to calculate, 

dF vac 1 

(53) 

8F 

v r vac A 

Q — A ex it 

vPamb 

(54) 

dF 

ujr vac 

0 4 — Pamb 

uA ex it 

(55) 


Eq. (53) reflects the fact that a unit change in site thrust is directly converted to a unit change 
in vacuum thrust. However, Eqs. (54) shows that a 1 psia change in ambient pressure causes a 
816 Ibf change in vacuum thrust. Likewise, a 1 inch 2 change in nozzle exit area causes a 0.036 Ibf 
change in vacuum thrust. Since the quantities used to compute the vacuum thrust were all measured 
quantities, there was no reason to calculate more partials derivatives using the chain rule as was the 
case of the venturi throat velocity and the mass flow quantities. 

Vacuum specific impulse. Two influence coefficients result from the specific impulse calculation. 


di 

dF va 


*p _ 


8L 


*T _ 


171TOT 

F vac 


draper 


. 2 

m TOT 


(56) 


(57) 


Typical values are 0.90 f fj and -0.43 lh 'j C sec for Eqs. (56) and (57), respectively. To express the 
partials with respect to measured quantities, the following equations were used: 


dl tp __ dl $p draper 

dpitnc draper &PUnt 

dlsp _ drh TOT 

dTu ne draper dTiint 


(58) 

(59) 
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81 ., = 

81 ., 

drriTOT 

(60) 

ddp 

dtriToT 

ddp 

dl.p 

dl. p 

dmj>OT 

(61) 

ddlint 

dm-roT 

ddline 

II 

dl.p 

drriTOT 

(62) 

adth 

dmTQT 

ddth 

dl.p 

II 

dF vac 

(63) 

dF, lU 

&Fvac 

dF.ite 

81 ., 

II 

> 

ys 

dF vac 

(64) 

dPamb 

dFvac 

dpamb 

dl, p 

II 

! Q) 

dF vac 

(65) 


dF vac 

dA txxt 


Values for the line pressure partial derivative were -0.3 and -0.1 for the oxygen and fuel 
sides, respectively. Oxygen and hydrogen line temperature partials took on values of 0.40 and 0.12 
Specific impulse partial derivatives with respect to differential pressure had values of -6.66 and 
-1.39 for the oxidizer and fuel sides, respectively with units of ~*f a . The partial with respect to line 
diameter had values of 1.67 and 0.78 and for the venturi throat diameter, -1732.0 and -577.9 

» ^ or ^e oxygen and hydrogen lines, respectively. A unit change in the site thrust value had a 
0.9 effect on specific impulse. The partial with respect to ambient pressure was 732.9 and 
with respect to nozzle exit area was 0.0325 . 

A sample calculation of the determination of the influence coefficient, error propagation, and 
final uncertainty computation for a simple case, the oxygen line contraction ratio, (3 ox , is given in 
Appendix B. 

Data Used in Study 

The data which were selected results for this type of experimental set-up was collected at the 
NASA RETF during the winter of 1986. Gaseous hydrogen and gaseous oxygen were the rocket 
engine’s respective fuel and oxidizer. Attached to the rocket engine combustion chamber was an 
optimized Rao contour nozzle which had been designed with an inviscid core area ratio of 1000:1 
and had been boundary layer corrected to an area ratio of 1030:1. 

Elemental Error Sources 

Before propagating errors from different instruments through the data reduction equations using 
the appropriate influence coefficients, an exhaustive list of possible error sources was compiled. The 
list is categorized by measurement device (pressure transducer, temperature thermocouple, and force 
transducer) and is subdivided into calibration and data acquisition errors. 

Pressure measurements (Table IV). The elemental sources of error for pressure and differential 
pressure transducers are combined. Zero reading errors in all pressure transducers were eliminated 
by numerically zeroing the instruments before taking data. Traceability of the transducer calibration 
from the standards lab to the National Bureau of Standards (NBS) was accounted for with NBS 
calibration reports on lab standard instruments. Errors due to changes in transducer calibration 
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pressure were eliminated by calibrating the instruments at test conditions. Transducer hysteresis 
and nonlinearity were determined by independent calibration of the instruments. 

Possible sources of data acquisition errors in pressure measurements, aside from data scatter, were 
manifested in the following ways. Transducer temperature differences at the time of zero balancing 
and data-taking are considered negligible because the propellant lines were at room temperature at 
all times. Errors due to the determination of a reference pressure were determined by an independent 
calibration. Changes in line temperature and pressure were determined to have a negligible effect on 
the transducer because they were at room temperature and experienced very small pressure changes 
during the tests. Damping of the propellant lines eliminated the effect of vibration on the transducer. 

Temperature measurements (Table V). Calibration errors for temperature thermocouples were 
characterized by the following errors. Specifications by the manufacturer identified wire calibration 
and reference temperature level determination errors. 

Data acquisition errors of the temperature measurement system included the following elements. 
Errors due to thermocouple temperature difference were assumed to be negligible because the pro- 
pellant lines were at room temperature at all times. Fabrication of the thermocouples were according 
to accepted standard practices and it was assumed that no errors were introduced as a result. Ef- 
fects of vibration on the thermocouple were assumed to be negligible because the propellant lines 
were damped. Because line pressure changes were negligible during testing, errors due to this source 
were neglected. Reference temperature stability was specified by the manufacturer to be ± 0.25%. 
Errors in thermocouple design due to radiation, friction, etc., when measuring gas temperatures 
were assumed to be minimized due to good thermocouple design practices. Heat conduction error 
sources and errors due to temperature gradients along non-homogeneous thermocouple wire were 
assumed to be negligible. 

Force measurements (Table VI). Calibration error in force measurements, including the following 
sources, were assessed. Errors due to the standards lab calibration of the transducer, including NBS 
traceability, were determined to be the error of the standards lab’s calibration instruments. Thrust 
stand hysteresis and non-linearity were determined experimentally to be ± 0.5% full scale. Errors 
due to a shift in load cell calibration caused by the attachment of adaptors and flexures were 
compensated for by in situ calibrations. 

Data acquisition errors which were addressed are listed below. Zero reading error was eliminated 
because the instruments were numerically zeroed before data taking. Vibration effects on the load 
cell and on the thrust stand were detemined to be negligible due to the design of the facility. Errors 
resulting from the misalignment between the engine force vector and the force vector measured by 
the data load cell train were minimized by using a triad of transducers in a parallel configuration. 
Similarly, the triad arrangement of force transducers accounts for misalignment of forces on an axis 
different from the engine centerline. Errors introduced due to pressurization on the load cell were 
corrected by aneroid calibration (ref. 28). Cell pressure effects on the cell wall were neglected 
because the cell wall was not the ground for the load cell. Errors due to the effect of changes in line 
pressure on the tare forces exerted on the thrust measurement system by propellant lines routed to 
the engine are reduced to a negligible magnitude by the physical set-up of the experiment. Effect of 
changes in temperature on the load cell was negligible because temperature changes were very small 
during calibration and testing. Similarly, thermal growth of the thrust stand errors were negligible 
also. Errors due to secondary airflow effects on the load cell were negligible during steady state 
operation of the engine. 

Error sources com mon to all instruments. In common to all the measurement devices were the 
following error source considerations. Error in the ability to determine a representative data value 
over a specified time interval as the data varied was reduced by averaging over five samples taken at 
50 Hertz intervals and characterizing these averages by a line of best fit. Errors introduced through 
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signal conditioning, electrical calibration, and digital systems were approximated with data from 
other test facilities at NASA LeRC. 


RESULTS 

By applying the methodology described and through the use of the equations listed, an uncer- 
tainty analysis was performed for the rocket engine test program conducted at the NASA Lewis 
RETF. A listing of the parameters measured and calculated, their nominal value, the amount of 
precision error in each, their associated degrees of freedom, the value of uncertainty, and their uncer- 
tainty expressed as a percentage, are tabulated in table VII. Most worthy of notice is the uncertainty 
value of specific impulse as a percentage. The test facility configured as described in this report is 
capable of measuring the specific impulse within 1.30%. Other parameters, such as vacuum thrust 
and total mass flow rate, were measured to within 1.12% and 0.72%, respectively. 

In order to appreciate which error sources most contribute to the specific impulse error, each 
term of the overall specific impulse precision error equation is expanded in table VIII. Each mea- 
sured parameter’s calibration and data acquisition error value is listed (in columns 2 and 4) as well 
as the influence coefficient of specific impulse taken with respect to that measured parameter (col- 
umn 1). The product of the influence coefficient and the associated error component (calibration 
or data acquisition) results in the amount of error, in units of specific impulse, which can be at- 
tributed to each measured parameter (columns 3 and 5). In presenting the data in this manner, the 
largest contributors to overall specific impulse error are easily seen. Table VIII shows readily that 
the calibration error is the major portion of the overall specific impulse error, indicating that the 
calibration process should be improved to effect the largest gains in specific impulse measurement 
accuracy. Within the category of calibration errors, the largest component of error is due to the 
ambient pressure transducer with over 43% of the total specific impulse calibration error. Another 
26% is contributed by the measured site thrust measurement. The next measurement which most 
affects the specific impulse calibration error is the oxygen venturi throat diameter with 14%. There- 
fore, the calibration errors of these three measurement devices are the cause of approximately 83% 
of the total specific calibration error. The six largest contributors to calibration error combine to 
comprise approximately 90% of the total error. Hence, to improve the overall uncertainty of specific 
impulse, better calibration of the these three instruments would effect the largest improvement. 

CONCLUSIONS 

An uncertainty analysis was conducted for the High Area Ratio Nozzle test program tests which 
were run during the winter of 1986 at the NASA Lewis Rocket Engine Test Facility (RETF). The 
standardized uncertainty methodology, as described by Abernethy (ref. 1), was followed. The data 
reduction equations were described and the expressions which define the influence coefficients were 
explained in detail. Elemental error sources of all types of measurements were discussed. 

A comparison of influence coefficients as determined by this analysis with similar coefficients, 
calculated with the data reduction program used in the experiment, was used to verify the uncertainty 
analysis model. In this way, the use of influence coefficients which were not taken with respect to 
the originally measured parameters were determined to be in error for use with this analysis. Their 
use leads to largely inflated results and caution should be taken to verify the influence coefficients 
values with independent computer programs whenever possible. 

Total uncertainty on specific impulse was determined to be 1.30%. Uncertainty on total mass 
flow rate was 0.72% and on vacuum thrust was 1.12%. 
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To continue the study of uncertainty analysis as it applies to rocket engine tests, the inclusion of 
estimates for the small bias errors of unknown magnitude for all instrumentation would provide better 
estimates of uncertainty in the performance parameters of interest to the rocket engine community. 


APPENDIX A 
List of Symbols 


A Area, inches 2 

B Bias error 

C Number of constants in a curve- 

fit 

C d Venturi discharge coefficient 

d Diameter, inches 

dp Static pressure difference be- 

tween the propellant line and the 
venturi throat 
F Thrust, Ibf 

h Enthalpy, BTU/lbm 

L P Specific impulse, sec 

K Conversion factor, 50079.5 

[(ft/secf/(BTU/lb m )} 

m Mass flow rate, Ibm/sec 

N Number of points in test 

P Parameters 

p Static pressure, psia 

r Computed result 

5 Precision error 


s Entropy, BTU/lb m /R 

SEE Standard estimate of error 

T Temperature, R 

<95 Student’s t value 

U 95 Uncertainty range encompassing 

95% of data 

U 99 Uncertainty range encompassing 

99% of data 

V Velocity, ft/ sec 

X Variable 

A' Average over all X values 

y Parameter 

/ 3 Propellant line contraction ratio 

9 Influence Coefficients 

v Degrees of freedom 

ir 3.14159 

p Density, lb m /ft 3 

<f> Generalized thermodynamic 

property variable 


Subscripts 


c 

Calculated 

site 

Site 

cal 

Calibration 

th 

Venturi throat station 

da 

Data Acquisition 

TOT 

Total (i.e M combined propellants) 

exit 

Nozzle exit plane station 

vac 

Vacuum 

i 

index 

P 

Pooled 

line 

Propellant line station 

X 

Variable 
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APPENDIX B 


Sample Calculation 

To demonstrate the calculation procedure, a simple example was chosen for which the given data 
is presented, the influence coefficients are calculated, the errors are propagated, and the uncertainty 
is determined. The oxygen line contraction ratio was chosen to illustrate the computation procedures 
discussed in this report. 

Given D ata: 

— 1.61 inches 
d t h y0 x = 0.406 inches 

v d, in . ... = 30 
= 30 

S ca \ dt - 0.0001 inches 

co*i«hne,or 

S C ai,d th tt T = 0.00065 inches 
s da,d linc „ T = 0 inches 
Sda,d lh „ r = 0 inches 

Calculation of Contraction Ratio, (3 


Pox 


dth,oa 

dhne.ox 


Pox 


(0.406) 

(1.61) 


0.252 


Calculation of Influence Coefficients 


dfiox _ d t h ,ox 

ddn neox df ine ox 


&0OX 

ddn neiOX 


(0.406) 

(1.61)2 


—0.157 inches' 1 


d/3 ox 1 

dd th ,ox dii ne ,ox 


dp ox 

ddiine t o x 


— - — rr 0.621 inches 1 
1.61 
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Propagation of Precision Errors 
Calibration Errors: 


Scal,0 ow 


d(3 a 


{ ddi 


line. ox 


-S ca i,d, in ...) 2 + ( 




dd 


th,oz 


■ S C al,d tk „ m ) 2 


S cattl ,, t = \/(— 0.157 x 0.0001) 3 + (0.621 x 0.00065) 2 
Scalp.. = v/(2.47x 10-10) + (1.63 x 10 -T) 


5e.w„ = 4.04 x 10" 4 


Data Acquisition Errors: 


S da,0„ 



line, ox 


S da,d lin + ( 


dPo: 


dd 


th,ox 


~ Sda,dth,*w ) 


S da< ^ m = v/(-0.157 x 0 + (0.621 x 0) 

Sda,0„ T = 0 


Combination of Errors: 


5 0». “ \/ 5 ?at ,/3„. + S da,0„ r 

Sp ot = v /( 4.04 x 10“ 4 ) 2 + 0 J 
S fi . m = 4.04 x IQ" 4 
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Calculation of the Degrees of Freedom 


(eL s?) : 


U 0om N S*_ 

2-/i=l Vi 


\212 


*'/9« = 




[(-0.157 x 0.0001) 2 + (0.621 x 0.00065)*] 

“ [(-0.T^7T0.^0T)5/W+ (0.621 x 0.00065) 4 /30] 


_ [2.47 x 1Q~ 10 + 1.63 x 10~ 7 ] 2 
"P"' ~ [2.03 x 10- 21 + 8.87 x 10" 16 ] 


2.665 x 10 


-14 


Vp ” 8.87 x 10- 16 


= 30.06 = 30 


Calculation of Uncertainty 

From Table I, the Student’s t value corresponding to vp or = 30 is < 95 = 2. 


Sp or — ±2 x (4.04 x 10 4 ) 


U 99,/3 0 * 


= 8.08 x 10“ 4 


To express the uncertainty as a percentage of the nominal value: 


^99 t p or % “ 


U M,p or 

Pox 


x 100 


U 99 } p„% 


8.08 x lOj^ 4 
0.252 


x 100 = 0.321% 
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TABLE I. DEGREES OF FREEDOM, v, AND ASSOCIATED 
STUDENT’S T VALUE, t 95 


H 

*95 

ID 

*95 

in 

*95 

ID 

*95 

ID 

*95 

V 

*95 

1 


mm 

2.447 

B 

EBTiTB 

16 

ESB1 

21 


26 

2.056 

E 

4.303 

B 

2.365 

m 

2.179 

17 

HBfl 

22 


27 


H 

3.182 

8 

2.306 

B 


18 

2.101 

23 


28 

2.048 

□ 

2.776 

m 

2.262 

14 

2.145 

19 


24 


29 

wm 

B 

2.571 

B 

2.228 

15 

2.131 

20 


25 


IBS 

mi 


TABLE n. INFLUENCE COEFFICIENTS OF 
THERMODYNAMIC PROPERTIES 
IN THE PROPELLANT LINE 


Influence Coefficients 

Oxygen Line 

Hydrogen Line 

Units 

^Itne /spline 

-0.790 x 10~ 2 

0.104 x 10- 1 

BTU/lb in 
tbf/in 2 — 

dhn ne / dTu ne 

0.239 x 10 +o 


BTU/lb,„ 

R 

dpiine i Spline 

0.622 x 10” J 

0.335 x 10“ 3 

HEnroiall 

dplinc / dT\ int 

-0.896 x 10~ 2 

-0.406 x lO' 3 

R 

dSlint / dpiine 

-0.107 x 10~ 3 

-0.156 x 10" 2 

ImSSill 

dsiine / dTline 

0.467 x lO" 3 

0.671 x 10" 2 

BTU/lbtnlR 
lb f /tn 3 


TABLE III. INFLUENCE COEFFICIENTS OF 
THERMODYNAMIC PROPERTIES 
AT THE VENTURI THROAT 


Influence Coefficients 

Oxygen Line 

Hydrogen Line 

Units 

dh th /dpth 

0.485 x 10' 1 

0.886 x 10 +o 

WFU/lb,„ 
lb ,fin * 

dhth/ dsii ne 

0.505 x 10+ 3 

0.524 x 10+ 3 

BTU/lb,,. 

BTUflb„JR 

dpth/dpth 



IK, /ft* 

UulifLl 

dpth/ dsn ne 

-0.187 x 10+ 2 


BTUJlb m J_R 

dTth/dpth 

0.238 x 10+° 


R 

lb f jin* 

dT th /da lint 

0.212 x 10+ 4 

0.146 x 10+ 3 

R 

BTU/lb„,/R 
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TABLE IV. SOURCES OF ELEMENTAL ERROR IN 
PRESSURE MEASUREMENTS 


Calibration errors due to... 

Comments 

Zero reading error 

Instruments are zeroed before data-taking 

Standards lab calibration of the trans- 
ducer calibration, including NBS trace- 
ability 

Errors of standards lab’s instruments 

Changes in transducer calibration pres- 
sure 

Not applicable 

Transducer hysteresis 

Independent calibration 

Transducer nonlinearity 

Independent calibration j 

Data acquisition errors due to... 

Comments 

Transducer temperature difference at zero 
balance 

Negligible because lines were at room 
temperature 

Transducer temperature difference at 
data-taking 

Negligible because lines were at room 
temperature 

Determination of reference pressure 

Independent calibration 

Changes in temperature of the transducer 

Negligible because lines were at room 
temperature 

Effect of vibration on the transducer 

Negligible because propellant lines were 
damped 

Effect of changes in line pressure 

Line pressure changes were negligible dur- 
ing testing 

Ability to determine a representative 
value over a specified time interval as data 
varies 

Data averaging reduces error to negligible 
magnitude 

Signal conditioning, electrical calibra- 
tions, and digital system 

Based on past data from other test facili- 
ties at NASA Lewis 
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TABLE V. SOURCES OF ELEMENTAL ERROR IN 
TEMPERATURE MEASUREMENTS 


Calibration errors due to... 

Comments 

Manufacturer’s specification of wire cali- 
bration 

Manufacturer’s specifications 

Reference temperature level detemination 

Manufacturer’s specifications 

Data acquisition errors due to... 

Comments 

T/C temperature difference 

Negligible because lines were at room 
temperature 

Fabrication of T/C 

As per standards 

Effect of vibration on the T/C 

Negligible because propellant lines are 
damped 

Effect of changes in line pressure 

Line pressure changes were during testing 
negligible 

Reference temperature stability 

±0.25% 

T/C design due to radiation, friction, etc., 
when measuring gas temperature 

Good design practice reduces this to neg- 
ligible magnitude 

Heat conduction 

All parts at uniform temperature 

Temperature gradients along nonhomoge- 
neous T/C wire 

All parts at uniform temperature 

Ability to determine a representative val- 
ues over a specified time interval as data 
varies 

Data averaging reduces error to negligible 
magnitude 

Signal conditioning, electrical calibra- 
tions, and digital system 

Based on past data from other test facili- 
ties at NASA Lewis 
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TABLE VI. SOURCES OF ELEMENTAL ERROR IN THRUST MEASUREMENTS 


Calibration errors due to... 

Comments 

Standards lab calibration of the transducer 
calibration, including traceability to NBS 

Errors of standards lab’s 
instruments 

Thrust stand hysteresis and non-linearity 

±0.5% as per manufacturer’s specifications 

Shift in load cell calibration caused by 
attachement of adaptors/flexures 

Compensated for with 
in situ calibrations 

Data acquisition errors due to... 

Comments 

Zero reading error 

Instruments are numerically 
zeroed before data-taking 

Effect of vibration on the load cell 

Negligible 

Effect of vibration on the thrust stand 

Negligible 

Misalignment between the engine force vector and 
the force vector measured by the load cell train 

1 

Triad arrangement accounts for 
misalignment 

Measurement of forces on an axis different from 
the engine centerline 

Triad arrangement accounts for misalignment 

Pressurization on the load cell 

Corrected by aneroid calibration 

Effect of changes in cell pressure on the test 
cell wall 

i 

Cell wall is not ground for load cell 

Effect of changes in line pressure on tare 
exerted on thrust measurement system by 
propellant lines, etc., routed to the engine 

Physical set-up reduces these 
to negligible magnitude 

Effect of changes in temperature on the 
load cell 

Temperature changes were negligible 
during calibration and testing 

Thermal growth of the thrust stand 

Negligible 

Secondary airflow effect on the load cell 

Negligible during steady-state operation 

Ability to determine a representative value 
over a specified time interval as data varies 

Data averaging reduces error 
to negligible magnitude 

Signal conditioning, electrical calibrations, 
and digital system 

Based on past data from other test 
facilities at NASA Lewis 
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TABLE VII. UNCERTAINTY OF MEASURED AND CALCULATED PARAMETERS 
(Part 1 of 2) 


Parameter 

Name 

[units] 

Nominal 

Value 

in[unit$] 

(i) 

Preci- 

sion 

Index (S) 
in [units] 
(2) 

Deg. 

of 

Free 

V 

(3) 

Uncrtnty 
Value 
<05 X (2) 
in [units] 

(4) 

Uncrtnty 
as % of 
Norn Val 

(i) * 100 
(5) 

Oxygen Line Pressure 

(pjia) 

649.40 

1.567 

5 

4.028 

0.620 

Hydrogen Line Pressure 

(psia) 

639.30 

7.275 

5 

18.704 

2.926 

Oxygen Line Temperature 

(*) 

510.80 

0.078 

5 

0.202 

0.039 

Hydrogen Line Temperature 

(*) 

533.60 

0.078 

5 

0.201 

0.038 

Oxygen Delta Pressure #1 

(pm) 

25.12 

0.136 

5 

0.350 

1.394 

Hydrogen Delta Pressure #1 

(psia) 

41.57 

0.302 

5 

0.775 

1.865 

Oxygen Delta Pressure #2 

(psia) 

25.26 

0.137 

5 

0.351 

1.389 

Hydrogen Delta Pressure #2 

(psia) 

42.16 

0.302 

5 

0.776 

1.839 

Oxygen Line Diameter 

(inches) 

1.61 

0.000 

30 

0.000 

0.000 

Hydrogen Line Diameter 

(inches) 

1.61 

0.000 

30 

0.000 

0.000 

Oxygen Venturi Throat Diameter 

(inches) 

0.41 

0.001 

30 

0.001 

0.320 

Hydrogen Venturi Throat Diameter 

(inches) 

0.44 

0.001 

30 

0.002 

0.343 

Site Thrust 

(»/) 

500.38 

1.684 

5 

4.329 

0.865 

Nozzle Exit Area 

(inches 2 ) 

815.70 

0.010 

30 

0.020 

0.002 

Ambient Pressure #1 

(psia) 

0.03 

0.004 

5 

0.010 

29.500 

Ambient Pressure #2 

(psia) 

0.04 

0.004 

5 

0.010 

24.378 

Average Ambient Pressure 

(psia) 

0.04 

0.004 

5 

0.010 

26.724 

Oxygen Line Contraction Ratio 

none 

0.25 

4 x 10- 4 

30 

0.001 

0.320 

Hydrogen Line Contraction Ratio 

none 

0.27 

5 x 10" 4 

30 

0.001 

0.343 

Oxygen Line Density 

(Ibm/ft 3 ) 

3.92 

0.010 

5 

0.026 

0.653 

Hydrogen Line Density 

(Ibjft 3 ) 

0.22 

0.002 

5 

0.006 

2.863 

Oxygen Line Entropy 

( BTU/lb m /R ) 

1.28 

2 x 10~ 4 

5 

0.001 

0.036 

Hydrogen Line Entropy 

(BTV/lb m / R) 

11.72 

0.011 

5 

0.029 

0.251 

Oxygen Line Enthalpy 

(BTU/lb m ) 

224.09 

0.026 

8 

0.061 

0.027 

Hydrogen Line Enthalpy 

(BTU/l b m ) 

1790.35 

0.349 

5 

0.898 

0.050 
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TABLE VII. UNCERTAINTY OF MEASURED AND CALCULATED PARAMETERS 
(Part 2 of 2) 


Parameter 

Name 

[units] 

Nominal 

Value 

in[units) 

(i) 

Preci- 

sion 

Index (S) 
in [units] 
(2) 

Deg. 

of 

Free 

V 

(3) 

Uncrtnty 

Value 

*95 X (2) 
in [units] 

(4) 

Uncrtnty 
as % of 
Norn Val 

(fj x ioo 

(5) 

Oxygen Throat Pressure #1 

(psia) 

624.28 

1.602 

5 

4.119 

0.660 

Hydrogen Throat Pressure #1 

(psia) 

597.73 

7.300 

5 

18.769 

3.140 

Oxygen Throat Density #1 

(«W ft 3 ) 

3.82 

0.010 

5 

0.026 

0.673 

Hydrogen Throat Density #1 

(«W ft 3 ) 

0.21 

0.002 

5 

0.006 

3.012 

Oxygen Throat Temperature #1 

(R) 

504.90 

0.104 

7 

0.245 

0.049 

Hydrogen Throat Temperature #1 

(R) 

523.71 

0.163 

12 

0.356 

0.068 

Oxygen Throat Enthalpy #1 

(BTU/lbm) 

222.88 

0.026 

8 

0.060 

0.027 

Hydrogen Throat Enthalpy #1 

(BTU/lb m ) 

1754.42 

0.642 

11 

1.413 

0.081 

Oxygen Throat Pressure #2 

(psia) 

624.14 

1.602 

5 

4.120 

0.660 

Hydrogen Throat Pressure #2 

(psia) 

597.14 

7.300 

5 

18.769 

3.143 

Oxygen Throat Density #2 

(«W ft 3 ) 

3.82 

0.010 

5 

0.026 

0.673 

Hydrogen Throat Density #2 

(lbm/ 'ft 3 ) 

0.21 

0.002 

5 

0.006 

3.014 

Oxygen Throat Temperature #2 

(R) 

504.87 

0.104 

7 

0.246 

0.049 

Hydrogen Throat Temperature #2 

(R) 

523.57 

0.164 

12 

0.358 

0.068 

Oxygen Throat Enthalpy #2 

(BTU/lb m ) 

222.88 

0.026 

8 

0.060 

0.027 

Hydrogen Throat Enthalpy #2 

(BTU/lb m ) 

1753.90 

0.646 

11 

1.421 

0.081 

Oxygen Throat Velocity #1 

(ft /sec) 

245.73 

0.850 

7 

2.009 

0.818 

Hydrogen Throat Velocity #1 

(ft/ sec) 

1344.67 

9.043 

8 

20.854 

1.551 

Oxygen Throat Velocity #2 

(ft/sec) 

246.43 

0.856 

6 

2.094 

0.850 

Hydrogen Throat Velocity #2 

(ft/ sec) 

1354.40 

9.055 

8 

20.880 

1.542 

Oxygen Mass Flow Rate #1 

(lb m / sec) 

0.82 

0.004 

26 

0.008 

0.949 

Hydrogen Mass Flow Rate #1 

(lbm /sec) 

0.29 

0.002 

n 

0.005 

1.728 

Oxygen Mass Flow Rate #2 

(lbm/ sec) 

0.82 

0.003 

26 

0.008 

0.953 

Hydrogen Mass Flow Rate #2 

(lbm /sec) 

0.30 

0.002 

11 

0.005 

1.728 

Oxygen Average Mass Flow Rate 

(lbm /sec) 

0.82 

0.003 

45 

0.007 

0.819 

Hydrogen Average Mass Flow Rate 

(lbm /sec) 

0.30 

0.002 

10 

0.005 

1.667 

Total Mass Flow Rate 

(lbm /sec) 

1.11 

0.004 

49 

0.008 

0.721 

Vacuum Thrust 

(Ibf) 

529.89 

2.771 

14 

5.943 

1.122 

Vacuum Specific Impulse 

(sec) 

476.10 

3.023 

29 

6.183 

1.299 
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TABLE VIII. SPECIFIC IMPULSE PRECISION ERROR BREAKDOWN 


Parameter 

Name 


Influence 
Coefficient 
sec/[units ] 

(i) 

Calib- 

ration 

Error 

[units] 

(2) 

Calib- 

ration 

Value 

[sec] 

1(1) x (2)1 
(3) 

Data 

Acqu. 

Error 

[units] 

(4) 

Data 

Acqu. 

Value 

[sec] 

1(1) x (4)| 
(5) 

Ambient Pressure #1 

(paia) 

366.40 

0.004 

1.376 

2 x 10~4 

0.064 

Ambient Pressure #2 

(psia) 

366.40 

0.004 

1.376 

2 X 10- 4 

0.064 

Site Thrust 

(»/) 

0.8985 

1.677 

1.498 

0.399 

0.358 

Uxygen Throat Diameter 

(inches) 

-1733.0 

7 X 1U 4 

T7TZ6 

0.000 

0.000 

Fuel Line Pressure 

(para) 

-0.1069 

7.272 

0.777 

0.566 

0.061 

Oxygen Line Pressure 

(psia) 

-0.0294 

1.562 

0.459 

0.318 

0.094 

Oxygen Delta Pressure #1 

(psia) 

-3.341 

0.132 

0.441 

0.088 

0.295 

Oxygen Delta Pressure #2 

(psia) 

-3.329 

0.132 

0.439 

0.092 

0.306 

Fuel Throat Diameter 

(inches) 

-577.90 

8 x 10 -4 

0.433 

0.000 

0.000 

Fuel Delta Pressure #1 

(psia) 

-0.697 

0.301 

0.210 

0.049 

0.034 

Fuel Delta Pressure #2 

(psia) 

-0.691 

0.301 

0.208 

0.051 

0.035 

Oxygen Line Temperature 

(R) 

0.3995 

0.075 

0.030 

0.061 

0.024 

Fuel Line Temperature 

(R) 

0.1160 

0.075 

0.009 

0.059 

0.007 

Nozzle Exit Area 

(inches 2 ) 

0.0325 

0.010 

3 x lO" 4 

6 x 10“ 5 

2 x 10~ 5 

Oxygen Line Diameter 

(inches) 

1.6730 

0.000 

0.000 

0.000 

0.000 

Fuel Line Diameter 

(inches) 

0.7787 

0.000 

0.000 

0.000 

0.000 

Root Sum Squared (RSS) 






±0.578 
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